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We study the dynamics of the universe with a scalar field and an SU(2) non-Abelian Gauge 
(Yang-Mills) field. The scalar field has an exponential potential and the Yang-Mills field is coupled 
to the scalar field with an exponential function of the scalar field. We find that the magnetic 
component of the Yang-Mills field assists acceleration of the cosmic expansion and a power-law 
inflation becomes possible even if the scalar field potential is steep, which may be expected from some 
compactification of higher-dimensional unified theories of fundamental interactions. This power-law 
inflationary solution is a stable attractor in a certain range of coupling parameters. Unlike the case 
with multiple Abelian gauge fields, the power-law inflationary solution with the dominant electric 
component is unstable because of the existence of non-linear coupling of the Yang-Mills field. We 
also analyze the dynamics for the non-inflationary regime, and find several attractor solutions. 



I. INTRODUCTION 



The idea of inflation now gives a standard sce- 
nario of the early evolution of the universe p]-@. 
It solves several difficulties such as the horizon 
and flatness problems in the Big-Bang cosmology, 
which has been confirmed by the precision cosmo- 
logical observations. It also provides us with a pre- 
diction on the origin of the observed density fluc- 
tuations. Many cosmological models with such a 
phase of accelerated expansion have been proposed 
by introducing a scalar field with an appropriate 
potential (or some alternative fields). However, it 
is desirable to derive a natural model from a fun- 
damental theory of particle physics without intro- 
ducing any anonymous field by hand. The most 
promising candidate for such a fundamental the- 
ory is the ten-dimensional superstring theory @ or 
eleven-dimensional M-theory@. They are hoped 
to give an interesting explanation for the acceler- 
ated expansion of the universe upon compactifica- 
tion to four dimensions. 

In the low-energy effective field theories of 
supcrstrings or supergravity theories, however, 
there is the so-called no-go theorem, which 
forbids such an inflating solution if the internal 
space is a time-independent nonsingular compact 
manifold without boundary Q. In order to evade 
this theorem, we have to violate some of those 
assumptions. We have three possibilities: 



'Electronic address: maeda" at" wascd a.jp| 

^Electronic address: K.Yamamoto"at"damtp.ca m.ac.uk| 



• a time-dependent internal space such as S- 
branes[9r{i2 

• an introduction of "singularity" such as 
branes [H0 

• a modification of gravitational action such as 
higher-curvature terms 0, [l6l - |20| 

Although some models could be promising, 
many models are still suffering from instability of 
a dilaton field or moduli fields. In fact, we natu- 
rally expect exponential couplings of moduli fields. 
Without fixing those moduli, many inflationary 
models are spoiled. 

An exponential coupling is not always harmful 
for inflation, however. In fact, we can find a power- 
law inflationplj with an exponential potential (22l 
l23l ]. It also provides the cosmic no hair theo- 
rem similar to the slow-roll inflation 2A\ . In su- 
pergravity theories and superstring models, an 
effective exponential potential Vq exp[— a<f>] natu- 
rally appears [25l - [27l ]. However, their potential is 
usually so steep that the power exponent of the 
scale factor cannot be much larger than unity, 
which makes it difficult to construct an accept- 
able inflationary model of the universe. For ex- 
ample, we find a = y2 and s/6 for two scalar 
fields in N = 2, six-dimensional supergravity 
model with S^-compactificationj^ . and the same 
is true for two scalar fields in N = 1, ten- 
dimensional supergravity model with gaugino con- 
densation [27} . Townsend summarized the possible 
exponential potentials derived by the compactifi- 
cation of ten- or eleven-dimensional supergravity 
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theories [28| . From flux compactifications, one ex- 
pects a > -s/6, while we may find y/2 < a < V& 
by hyperbolic compactifications. Neither of them 
offers a flat enough potential for inflation. 

In the unified theories of fundamental interac- 
tions, there naturally exist gauge fields, which 
may be included in the original action such as 
the heterotic string theory or can be induced 
by Kaluza-Klcin compactification. In effective 
four-dimensional theories derived from higher- 
dimensional unified theories, we also expect those 
gauge fields coupled exponentially to moduli fields 
such as \ exp[\(j}]F 2 . Hull and Townsend discussed 
such a coupling for the case of U(l) gauge fields. 
They found that the possible values of the cou- 
pling in the four-dimensional effective action are 
A = 0, \/2/3, y/2, or y/E in the context of black 
holes in the type II string theory compactified on 
a six torus [23. In M-thcory (eleven-dimensional 
supcrgravity) with intersecting branes, the four- 
dimensional effective action also contains the same 
moduli couplings to U(l) multiplet [30j ] . 

If the strengths of the couplings between gauge 
fields and a scalar field are similar to that of the 
scalar self-coupling, the gauge fields may affect the 
dynamics of the scalar field. In fact there are sev- 
eral discussions about the dynamics of inflation, 
where supportive roles of gauge fields in realizing 
accelerated expansion have been observed j3ll - l42T | . 

The effect of the gauge-kinetic coupling on the 
inflationary dynamics was first discussed in the 
context of anisotropic mflationpjll]. assuming a 
U(l) gauge field coupled to an inflaton field. Since 
a single U(l) field cannot exist in Friedmann- 
Lemaitre- Robertson- Walker (FLRW) isotropic and 
homogeneous spacetimes, they discussed Bianchi 
spacctimcs as the cosmological model. They spec- 
ified the scalar potential to be quadratic and chose 
exp[c0 2 ] as the gauge kinetic coupling. They 
showed that an anisotropic inflationary era may 
arise as a transient attractor state while the scalar 
inflaton is slowly rolling. The anisotropy even- 
tually disappears as the scalar field oscillates to- 
wards the end of inflation. The observational 
relic of the anisotropic inflationary era was also 
discussed [3ll HH . 

While the chaotic inflation driven by the 
quadratic potential is phcnomenologically interest- 
ing as it automatically results in reheating, the 
form of the inflaton-gauge interaction discussed in 
[3l| may not naturally appear in the unified theo- 
ries. They also studied the case with an exponen- 
tial potential and a U(l) gauge field coupled expo- 
nentially to the scalar field, which suits the frame- 
work of the unified theories better. They found an 
exact anisotropic inflationary solution, which is an 



attractor independent of the initial conditions [33j. 
Since our present universe is almost isotropic, this 
model must be severely constrained. 

However, if there exist more than two gauge 
fields, we find an interesting scenario. Although it 
requires an artificial assumption that all the gauge 
fields couple to the inflaton through a common 
gaugc-kinctic function, one can obtain a totally ho- 
mogeneous and isotropic inflationary solution as an 
attractor [Hj]. Since the anisotropic inflation can 
be found as a transient attractor, we might have 
a chance to find distinct observational signatures. 
An important result is that an isotropic power- 
law inflationary solution appears as an attractor 
even for a steep exponential potential for the infla- 
ton, which is expected from the unified theories of 
fundamental interactions. While there are certain 
conditions to be satisfied by the gauge-kinetic cou- 
pling constant, they are not so strict as the usual 
slow-roll conditions and could fall within the reach 
of the supergravity theories. 

In the case of U(l) multiplet fields, we usually 
expect the different gaugc-kinctic coupling con- 
stants for different fields in the context of the 
unified theories. However, if wc consider a non- 
Abclian gauge field, it consists of "multiple" vector 
fields with a single common gauge- kinetic coupling 
constant. As a result, the discouraging feature of 
U(l) multiplet will disappear. The conventional 
chaotic inflationary model with a non-Abelian 
gauge field has been studied 35]. Motivated by its 
phenomcnological development and the aforemen- 
tioned features of high-energy physics, in this pa- 
per, we study SU(2) non-Abelian gauge field cou- 
pled exponentially to a scalar field with an expo- 
nential potential, in order to know whether the 
non-Abelian gauge field has the similar nice prop- 
erties as the U(l) multiplet case. 

We should note that there is also an approach 
different from the present gaugc-kinctic coupling 
model[4l|, [42|. They consider an axion field cou- 
pled to a non-Abelian gauge field, which is named 
chromo-natural inflation. It may give another 
interesting inflationary regime with non-Abelian 
gauge fields. 

In the following, we present the basic equations 
of our system, and obtain power-law solutions in §. 
Mil Wc find that a power-law inflationary solution 
is found only for the case of the magnetic com- 
ponent dominance in contrast to the U(l) triplet 
case, in which both inflationary solutions with elec- 
tric field and magnetic field arc possible. In §. Hvl 
we describe the solutions as fixed points of a dy- 
namical system and analyze their stabilities. In §. 
IVl wc perform numerical analysis for the range of 
coupling constants where fixed points do not exist. 
It also tells us how the attractor state is achieved 
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from generic initial data. Concluding remarks and discussions are given in §. |VT 



J 



II. BASIC EQUATIONS 

We use the unit of k 2 = 8nG = 1. The action wc discuss is 
S = J d 4 xv^g 
where 

p(a) _ ft 4(a) _a /l(a) 1 „ , A{b) a(c) 



\r i(V0) 2 - V(cf>) \f{<t>)F$F^ 



is an SU(2) non-Abclian gauge field, which we call the Yang-Mills (YM) field, and g YM is its coupling 
constant. The coupling to the scalar field /(</>) and th e scalar potential V{4>) are given respectively by 



V(<f>) 



V e 



a can be set non-negative without loss of general- 
ity. We also restrict ourselves to Vq > since our 
primary interest here is inflation. 

Throughout the article, we discuss a flat FLRW 
spacetime[43l]. whose metric is given by 



ds 1 



-dt 2 + a 2 (t)d 2 x. 



We assume that the vector potential is given by 



-4 



(a) 



0,A 



(a) _ 



A{t)S. 



(a) 



(2.1) 



so that the YM field is taken to be isotropic. This 
configulation results in both homogeneous electric 
and magnetic components, which are written in 
the coordinate basis as 



E 



(a) 



F^=a(t)E(t)S^ 



*(*) 



with 



.4 
(1 



E := , and B = c/ym — ~ 



A 2 _ 

2 



being their comoving field strengths. This is an im- 
portant difference from U(l) gauge fields, for which 
we find only the electric component in the above 
vector potential. The homogeneous magnetic com- 
ponent in U(l) gauge fields is obtained only when 
we introduce an appropriate inhomogeneous vector 
potential. As a result, the electric component and 
magnetic one in the U(l) fields are independent. 
We can discuss each component separately. In con- 
trast, the YM field always consists of two compo- 
nents in the above isotropic configuration (12.11) and 



the homogeneous field is found only by a homoge- 
neous vector potential. If one introduces any spa- 
tial dependence to the vector potential, the field 
strengths become inhomogeneous. We should also 
note that we need more than two U(l) fields with 
a common coupling to the scalar field as discussed 
in [HI in order to find an isotropic and homoge- 
neous attractor spacetime. Otherwise, we find an 
anisotropic universe. For an SU(2) gauge field, this 
uniform coupling is a necessary consequence of the 
symmetry. 

The Einstein equations are 
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H = - 



l -4> 2 + v + P , 



1 



(2.2) 
(2.3) 



where the dots denote the time derivative d/dt, and 
H = a/a is the Hubble expansion rate. p Ym is the 
YM energy density, which consists of the electric 
and magnetic parts, i.e., 



Pym = Pe+ Pb ■ 
They arc defined by 



PE 



-e^E 2 , p B = °-e x *B 2 . 



(2.4) 



(2.5) 



The equation of motion for the scalar field is 
'<j> + mj)-aV - \{pe-Pb) =0, (2.6) 

and the equation of motion for the YM field is 
simply 



A 3 

A + HA + \(j)A + 2g YM ^ = . 



(2.7) 



Using the Bianchi identity, Eq. (|2.3|) is obtained 
from Eqs. (J23), (HU) and (gT]). Hence wc take 
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and (|2.7[) as the basic equations of our 
system. 

The YM equation can be reduced to the first 
order equations for each energy density as 



A 

p E = -(4ff + \<t>)p E -A—pb 



A 

p B = ~(4ff - X(j))pB + 4-p B 



(21 



The terms with AH come from the radiation-like 
behavior of the YM field (p ra d oc a~ 4 ) and the 
last terms are from the non-linear interaction in 
the YM field. In fact, for U(l) triplet fields with 
a uniform exponential coupling to a scalar field, 
we find the evolution equations for energy densi- 
ties by dropping the non-linear interaction terms. 
As we shall see later, the non-linear terms play a 
significant role in the dynamics of the model. 



III. POWER-LAW SOLUTIONS 

Since we have the exponential potential, we ex- 
pect a power-law expansion and look for the possi- 
bility of power-law inflation. Suppose our solution 
is given by 



a t p 
2 



hit- 



(3.1) 
(3.2) 



where p is assumed to be a constant, and ao and 4>q 
are initial values. The coefficient 2 /a in front of In t 
is determined by requiring that the t dependence 
of <fi 2 and V be the same, in order to satisfy the 
Hamiltonian constraint (12.21). 



A. The case with U(l) triplet fields 



First 



consider the case 
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fields, which was discussed in 
to be solved are 



Pe = -(4:H + \<P)p E 
p B = -(4H - \(j>)p B : 



with U(l) triplet 
The equations 



(3.3) 



and Eqs. Q2JJ) and ([2T6]> . In our setting (|2Tj) . the 
magnetic field vanishes. However, if we add an ap- 
propriate inhomegeneoes vector potential, a home- 
geneous magnetic field can appear and the energy 
densities of the electromagnetic fields satisfy Eqs. 



Dynamics of the scalar field 



Eqs. (|3.3p are eaily integrated as 



Pe = Peo 



e -\(4>~4>o) 
(a/a ) 4 

e A(0-0o) 



PB = PBO- 



a/ao) 4 



(3.4) 



where PeoiPbot4>o & n d ao are integration con- 
stants. 

Then the equation of the scalar field is reduced 

to 



3*^ = 0, 



(3.5) 



where 

V eS := V e- a * + i (C E e-^ + C B e X4> ) (3.6) 
with 



a 1 



C E = PEoay* 



C 



B 



4 — A0o 

PBoa e 



(3.7) 



Although the original potential V is monotonically 
decreasing, the effective potential (|3.6[) has a mini- 
mum point for A < and Ce ^ 0, or for A > and 
Cb 0. As a result, the scalar field will evolve 
more slowly than the case only with the exponen- 
tial potential V. Since there exists the pre-factor 
a -4 , the minimum point will move and the min- 
imum value will decrease as the universe evolves. 
Hence we do not have an exponential expansion, 
but have a power-law expansion whose power ex- 
ponent is larger than the original power-law expan- 
sion driven solely by the potential V. This is the 
mechanism that a gauge field coupled to a scalar 
field assists slowing down the motion of the scalar 
field and inflationary expansion becomes possible 
even for a steep potential. 

Next we present the explicit power-law solutions. 
Assuming that the expansion of the universe and 
the evolution of the scalar field are described by 
Eqs. (|3.ip and p.2[) , and the energy densities of 
the electromagnetic fields are proportional to t~ 2 , 
i.e., pe = PEo/t 2 and ps = PBo/t 2 , we find that 
Eqs. (|3.3[) become the algebraic equations: 



2p + — ) 



Pbo 



2p ] Pbo 

a 



There are two cases :pso — an d Peo 
we shall discuss separately. 



(3-8) 
(3.9) 
0, which 
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2. The case with the electric field (-Bui) 

For the case with the dominant electric field 
(Pb = 0), which we shall call the regime -Bui) 
assuming pe = Peo^ 2 (peo : constant), we find 
three algebraic equations: Eq. (|3.8p and 

p2 = \(h + Vae ~ a *° + PE0 ) (3 - 10) 

1 aV 

a a 

which are rearranged into 

V e- a ^ = —r [4 - 3A(a - A)] (3.13) 
4ar 



Peq 



3 

4a 2 



a (a - A) - 4] 



(3-14) 



Since the left-hand-sides are positive definite, in 
order for such a solution to exist, we have to impose 
the following conditions: 



A < 



a - | {a < y/6) 

\(a- V« 2 -16/3) (a>V6) (3 ' 15) 



The power-law inflation is possible for the range 
of coupling parameters of A < —a and a (a — A) > 
4, as was shown in 13411 . 



3. The case with the magnetic field (-Bui) 

For the case only with the magnetic field {pE = 
0), which we shall call the regime -Bui, we find the 
same result by changing the sign of A. The solution 
is described by 



1 / A 
P=-(l + - 



^[4 + 3A(a + A)] (3.17) 



(3.16) 



Vne- 



3 



(3.18) 



and the existence conditions are 



A> 



-a + | (a < V6) 

- ^a 2 - 16/3) (a > V6) 



(3.19) 



The power-law inflation is obtained for the param- 
eter range of A > a and a(a + A) > 4. 

Defining the density parameters of each compo- 
nent by Qe = Pe/3H 2 (the electric field), = 
Pb/ZH 2 (the magnetic field), n v = V/3H 2 (the 



potential), and O^- = <f> 2 /6H 2 (the kinetic term of 
the scalar field), we find that those depend on the 
coupling parameters. We show one example for 
the power-law inflation with magnetic field in Fig. 
[TJ We find that the magnetic field gives a certain 
contribution to the expansion of the universe. 



(3.10) 


1.: 

08 




(3.11) 
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0.2 


- n K 


(3.12) 
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FIG. 1: The density parameters of each component 
[Qb =(the magentic field), ilv (the potential), and 
Q.K (the kinetic term of the scalar field)] for the case of 
a = v6 and A > a. The power exponent of the scale 
factor is given by p — (1 + A/a)/2. 



We will show the stability condition later in 

ttvU 



B. The case with YM field 

Now we show the dynamics changes when the 
YM interaction is turned on. Since we have the 
non-linear coupling in the YM field, a simple 
power-law ansatz may not work. But we first look 
for a solution similar to those found in the U(l) 
triplet case. 



1. The case with the dominant electric component 
(Bym) 

If we assume pe 3> Pb but ps being non- 
vanishing due to the interaction between electric 
and magnetic fields, which we shall call the regime 
E YM , dropping the term with pb, we find the same 
equations as the U(l) triplet case. As a result, we 
find the same solution [Eqs. (|3.12[) - (|3.14[) ] as long 
as the electric component stays dominant. Under 
the conditions A < —a and a(a— A) > 4, we obtain 
the power-law inflationary solution. Note that an 
accelerated expansion is possible even if a > \/2 
just as was the case for the U(l) triplet electric 
type inflation (33| . 

However, the situation in the case of YM field is 
not exactly the same as that for the U(l) triplet 
fields. In the above analysis, we have ignored the 
magnetic component, which is valid in the U(l) 
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triplet case because the electric and magnetic fields 
are decoupled. However, the electric and magnetic 
components are always coupled in the YM field. 
Then we have to check whether the magnetic com- 
ponent is always negligible or not when it is ini- 
tially small. 

Since we assume the magnetic component is ini- 
tially very small, we can solve the YM equation 
(|2.7|) dropping the term with g YM (the magnetic 
contribution) as 



A = Aia e 



A 



1 -A</> 



At a + 3A 
-t 5S~ 



a e 



Xtpo 



2aA 1 



(a - 3A)a e A *o 



where Aq and A\ are integration constants. Using 
this solution, we evaluate the ratio of two energy 
densities as 



PB_ 

PE 



Pb 

PE 



where pb/pe\o is the initial value. We drop Aq 
since we are interested in the asymptotic behavior 
(t — > oo). As a result, if the magnetic component 
is initially sufficiently small, we have a power-law 
inflation just as the case with the Abelian mul- 
tiple fields, but the contribution of the magnetic 
component gets larger during the evolution of the 
universe, and then the inflationary phase eventu- 
ally ends because of the growth of the magnetic 
component. The e- folding time when the approxi- 
mation becomes no longer valid is evaluated as 



-^-folding = ln(a e nd/«o) 
4 \PE 



For example, if we assume pb/pe\o = 10 -8 , we 
find JVe-foiding ~ 4.6. Hence unless the initial value 
of the magnetic energy density is extremely small, 
we do not find a sufficient e-folding number for the 
inflationary universe. 



2. The case with the dominant magnetic component 
(Bym) 

We shall call it the regime B YM , when the mag- 
netic component is much larger than the electric 
one. In this case the situation is not so simple as 
the U(l) case because we cannot ignore the non- 
linear term in (|2.8[) even if we assume pe <C Pb- 
Suppose we have the same solution as the U(l) 
case. We then evaluate A/ A by the YM equation 



(|2.7p . which is now 
A 



A + 2gl 



A'' 



2at 



We then find the asymptotic solution as 



A — Aqq 
which leads to 

A 

A 



Aa 2 a 2 A 

-tlx y YM ^-L 



2 

YM* 1 a OO 



t- 



(A — a)(a + 3X) a 2 



Aag 2 M A c 



(a 



3X)a 2 



-A/a 



(3.20) 



(3.21) 



(3.22) 



as t — > oo. This ratio decays faster than t 1 , at 
which rate AH ± X(f> evolve in Eq. (|2.8jl . So we can 
ignore the non- linear term with A/ A, which gives 
exactly the same equations as the U(l) magnetic 
case. 

As a result, we have the power-law solution just 
the same as in Bui- This solution is obtained 
asymptotically, and the non-linear term does not 
destroy it unlike the regime E YM where the electric 
components dominate. 



3. The case with both components 

If both electric component and magnetic one are 
of equal magnitude, we cannot ignore either of 
them. Since it is a non-linearly coupled system, 
it is difficult to figure out what kind of solutions 
we expect. Although we need numerical studies, 
which we will give later, here we shall discuss one 
simple case, in which we assume the power-law be- 
havior. 

Suppose that the YM potential A is a power- 
law function as A oc t s , and the scale factor and 
the scalar field are given by Eqs. (|3.ip and (j3.2|) . 
If the energy densities of the electric and magnetic 
fields arc similar, i.e., pe ~ Pb, we find s = p — 1, 
and then 



p E ~ Pb oc 



t -2 x t -2(i-A/«) 



Inserting this behavior into the YM equation (j2.7|) . 
we find 



2(p-l) p-1 



2gt 



= 0. 



which implies 



2(p-l) [P- 1 + -J <0- 

For the power-law inflation (p > 1), we have 
A 

1 > p > 1 . 

a 



7 



As a result, pe ~ Pb drops faster than cx t~ 2 , 
which is the scaling of energy density of the scalar 
field. Hence the contribution of YM field becomes 
less important as t — > oo. It appears that it is 
not possible to find a power-law inflation with a 
significant residual pe ~ Pb- Here, we find the 
power-law expansion only by a scalar field. An 
accelerated expansion is possible if a < y/2, as the 
conventional power-law inflation. 

In the next section, we will give more details of 
interesting solutions in the present system includ- 
ing the inflationary solutions we have found and 
analyze their stability as fixed points in a dynam- 
ical system. 



IV. STABILITY ANALYSIS 
A. Dynamical System 

In order to analyze the dynamical behavior of 
our solutions found in §. IIII B[ we rewrite the basic 
equations in the form of a first order autonomous 
system. The inflationary solutions discussed in 
the previous sections, along with other interest- 
ing ones, appear as fixed points in the dynamical 
system. This allows us to study their local sta- 
bility and reveal a complicated dynamical behav- 
ior that goes beyond the simple power-law time- 
dependence. We shall change the time coordinate 
from t to the e- folding number N = ln(a/ao), and 
introduce new variables normalized by the Hubble 
expansion rate H as 



£ 
B 
A 



A<t>_ = - e i\<t> 



E 

H 

D 



a 



pip— - A 2 



1/2 A, 

5ym e 4 



.4 



(4.1) 
(4.2) 
(4.3) 



Primes denote differentiations with respect to the 
e- folding number N. We then introduce the den- 
sity parameter of the YM field as 



3H 2 ~ 2 



(£ 2 + B 2 



(4.4) 



and those of the potential and the kinetic energy 
of the scalar field as 



n K = 



v 

'3H 2 ' 

2 _ vo 2 

m 2 ~ 6 



(4.5) 
(4.6) 



where w := (fi/H = ft. We also use 



A = 



PB - PE 

P\ r M 



(4.7) 



which describes the difference of the fractions of 
the magnetic and electric components. It en- 
ables a unified treatment of electric- and magnetic- 
dominant regimes and also makes the asymme- 
try clear when the YM coupling comes into play. 
A = 1 and —1 correspond to the regime B YM and 
E YM , respectively. 

The Friedmann equation (|2.2[) now reads 



V 



1 . 



(4.8) 



The equation for the scalar field (|2.6[) is 



w' = - (6 - vo 2 ) (a - w) + [2m - 3 (a + AA)] il YM . (4.9) 

where we have used the Friedmann equation (|4.8p to eliminate Vly. The equations for the YM field (|2.8[) 
are now 

A' = - Hot + A) - 4 (1 - fi YM )] A - T£ (4.10) 

£' = i[n7(tu-A)-4(l-f7 YM )]f + 2r^ 3 . (4.11) 

where 

T = g^e-^H-y\ (4.12) 
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whose evolution equation is 

T' = - [w{vj - A) + 4fi YM ] T . (4.13) 

This auxiliary quantity V is the "normalized" YM coupling in the sense that the subsystem defined by 
r = corresponds to the dynamical system that describes homogeneous and isotropic U(l) triplet fields. 
The YM equations (|4. 10[) and (|4. 1 1[) are rewritten in terms of the denisty parameter f2 YM and the ratio 
A as 

!1™ = [-4 + t t 7( G 7 + AA)+4f2 YM ]ft YM (4.14) 
A' = Aro(l - A 2 ) -4er(l - A)* (1 + A)* o| M , (4.15) 

where e = sign(.4£). 



We then find the dynamical system in a closed 
form. Since the physical interpretation of the nor- 
malized vector potential A is not clear, we take 
Eqs. gU), (|4T4) . (|4T5)) and (|443)) with the 
Hamiltonian constraint (= the Fricdmann equa- 
tion) (|4.8|) as the basic equations to analyze the 
stability around the fixed points. The drawback is 
the appearance of e which takes into account the 
ambiguity inherent to taking square roots. This 
causes a problem in the numerical study in the 
next section when the system undergoes oscilla- 
tions. For this reason, Eqs. (|4.10[) and (|4 . in- 
stead of Eqs. (|4.14[) and (|4.15|) are used there. 



B. The case with U(l) triplet fields 

Before going into analysis of our system, for an 
introduction and a comparison, we first summa- 
rize the case with the U(l) triplet fields, which 
was discussed in [34| . using the present dynamical 
variables. To make a clear distinction, we replace 
J7ym with fiui- Now A = 1 and —1 correspond to 
the regimes Btji & n d Em respectively. 

In the case with the U(l) triplet fields, the dy- 
namical system is obtained by setting T = in the 
above; 



1 



6 — n7 2 ) (a — w) 



O' 
A' 



+ [2m - 3 {a + AA)] n m 
[-4 + vo 2 + XrnA + 4f2ui] fl m 
Xvu (1 - A 2 ) . 



If zu ^ and Slui 7^ 0, the fixed points are 
classified into two cases; A = — 1 (the case with 
the electric field) and A = 1 (the case with the 
magnetic field). In each case, we find two fixed 
points as follows: 



(a) zu = — 3AA , f2ui 



2-3A 2 



(6) w 



AA 



f2ui = 



a(a + AA) - 4 
(a + AA) 2 



Since the density parameters arc positive defi- 
nite, A < \j2j'i for the fixed point (a) to exist. 
In this case, fly = 0, which means that either 
the potential is absent from the beginning or the 
potential becomes asymptotically negligible com- 
pared with the kinetic term vu 2 /2. 

From a perturbative analysis, we can check the 
stability of these fixed points. For the fixed points 
(a), we find that at least the eigenvalue for the 
perturbation of A is always positive (wa = 6A 2 ). 
Hence it is unstable. Hereafter, we use to to denote 
eigenvalues with subscripts indicating the variable 
to which the eigenvalue is associated. 

The fixed points (6) represent the power-law so- 
lutions {E\ji and Bui) found in §. IIII Al The per- 
turbative analysis gives the following three eigen- 
values: 



8AA 



a + AA 

and the two roots of the quadratic equation 



(4.16) 



(a + XAfiu 2 + (a + XA)(a + 3AA)w 
+[a(a + AA) - 4] [4 + 3AA(a + AA)] = , 

from the perturbations of w and f2 YM - 

From the existence conditions given by Eq. 
(pTT5)) or Eq. ([3~T^l) . we have a(a + AA) - 4 > 0. 
Hence, if and only if AA > is satisfied, all the 
three eigenvalues are negative (or the real parts 
are negative if they are complex). As a result, the 
solution with the conditions 



AA>0, a(a + AA)-4>0 



(4.17) 
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is stable against linear perturbations. More con- 
cretely, the stability conditions for the cases with 
the electric field (-Eui) and magnetic field (Bm) 
are given by, 

A < 0, a(a- A) > 4, (4.18) 
A > 0, a(a + A) > 4, (4.19) 

respectively. If A > 0, the magnetic power-law so- 
lution (Eqs. (J3IU, (GT2) and §^ with $Hty and 
(13 . 1 8[) ) is always an attractor in the parameter 
range of a(a + A) > 4, while if A < the elec- 
tric power-law solution (Eqs. (f3Tj) . (f3~2]) and f|3T4|) 
with p,12p and (|3.14[) ) is always an attractor in 
the parameter range of a (a — A) > 4. 

For the rest of the parameter space (a — 4/a < 
A < — ol + 4/a), the attractor is a fixed point with 
Qui = 0, where the scalar field dominates the uni- 
verse. The fixed point, which we denote 5tji, is 
given by 



The perturbative analysis gives the three eigenval- 



ues as 



w !2ui 



c?). 



a(a + AA) - 4, 
-2aXA. 



w = a, Qui =0, A = ±1 . 



(4.20) 



Hence the power-law solution driven only by the 
scalar field is stable if A A > and a(a + A A) — 
4 < 0. Between the two solutions with alternative 
sings, the stable one is 

zu = a, Qui=0, A = l. (4.21) 

for A > 0, while 

w = a, Oui=0, A = -l. (4.22) 

for A < 0. 



We summarize the result for the U(l) triplet case in Fig. [2] and in Table Q] 




FIG. 2: The parameter range for power-law solutions in the case with the U(l) triplet fields. The inflationary 
attractors are indicated by black letters while non-inflationary attractors are red. The attractor solution with 
the electric field is given for A < and a(a — A) > 4 (Em), while one with the magnetic field is for A > and 
a(a + A) > 4 (-Bui) ■ In the range of a — 4/a < A < —a + 4/a (Sui), we find the attractor dominated by the 
scalar field. All inflationary solutions are stable. The inflation with the U(l) triplet field is found in the range 
either of A < —a and a(a — A) > 4 (i?ui-I:with the electric field) or of A > a and a(a + A) > 4 (Bui-hwith the 
magnetic field). The conventional power law inflation with an exponential potential is possible only for a < y2 
and a — 4/a < A < -a + 4/a. (Sui-I) 



C. Important Fixed Points in the dynamical 
system with YM field 

Now we move on to include the non-linear YM 
interaction. The non-trivial fixed points are clas- 
sified into two cases: T = and r / 0. In the 
former case, we find the same fixed points as the 
U(f ) triplet case, although their stability is differ- 



ent as we will show later. The latter case gives new 
fixed points, which do not exist in the U(l) triplet 
system. 

Note that a fixed point may not be found by 
an exact solution, but can be reached as a certain 
limit. For example, the fixed points with T = 
would imply either <?ym = or He Xr ^/ 2 — oo, nei- 
ther of which is of interest in our analysis here. 
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fixed point 


Sui 


Evi 


Bui 


existence 


a < VE 


a — 4/a > A 


— a + 4/a < A 


P 


2/a 2 


| (1 - A/a) 


+ A/a) 







a (a — A) —4 
(a — A) 2 


a(a + A)— 4 
(a + Ap 




1 - a 2 /6 


4-3A(a-A) 
3(a-A)^ 


4+3A(a + A) 
3(a+A)^ 


Ok 


a 2 /6 


8 

3{a-\) A 


8 

3(a+A)^ 


stability 


a — 4/a < A < —a + 4/a 


A < 


A > 


inflation 


a < V2 (Sui-I) 


A < — a (£ui-I) 


A > a (Bm-I) 



TABLE I: The fixed points and their properties for the case with U(l) triplet fields. The second row gives the 
existence conditions. The bottom two rows are understood to hold when those existence conditions are satisfied. 



However, starting from gYM and finite H and 
<f>, the system may approach T — > asymptotically 
as t — > oo. From the mathematical point of view, 
those fixed points are well-defined and a part of 
the dynamical system and we include them in the 
following analysis. 

1. r = 

In this case, which should reproduce the fixed 
points of the previous subsection, we can classify 
the solutions into two cases: f2 Y M = and 

(a) » YM = 

In the case with S1 YM = 0, the scalar field energy 
is dominant. From (|4.9j) . we find either w 2 — 6 or 
w = a with A 2 = 1. The former fixed point cor- 
responds to the case that the kinetic energy of the 
scalar field is dominant, which is unstable against 
perturbations. The latter fixed points denote the 
power-law expanding universe with an exponen- 
tial potential (the counterpart of Sui) and will be 
called Sym- The ratio of the potential energy V to 
the kinetic energy is (6 — a 2 ) /a 2 . As is well known, 
these fixed points are attractors if a < y6 for the 
case only with a scalar field. 

In the present case, because of the YM field, 
the stability condition changes as follows. The lin- 
earized equations for these fixed points give four 
eigenvalues: 

^ = \ (« 2 - 6) , (4.23) 
^ = aAA + a 2 -4, (4.24) 
uja = -2aAA, (4.25) 

uj r = -a(a-A). (4.26) 

Three eigenvalues (gj25J) , P~2"I|) . and P~2"5j) are all 

negative if A > and A < —a + 4/a for A = 1, or 



if A < and A > a - 4/a for A = -1. The forth 
eigenvalue (|4.26|) becomes negative if A > a. As a 
result, the fixed point with A = 1 is stable in the 
parameter range ofa<A<— a + 4/a. 

On the other hand, taking A < a gives insta- 
bility against the perturbations of T. However, as 
long as Qym stays small, the growing T does not 
disturb the evolution of the universe as well as the 
dynamics of the scalar field since T docs not appear 
explicitly in Eqs. (|4.8[) and (|4.9j) . This is indeed 
the case when a — 4/a < A < a. As we shall con- 
firm later in the numerical analysis, the dynamics 
of the universe is dominated by the scalar field and 
accurately described by the fixed point discussed 
here, despite the apparent instability in the eigen- 
value wr ■ The exponentially increasing T only trig- 
gers a rapid oscillation for the perturbed YM field 
whose amplitude remains small. 

In summary, we conclude that these fixed points 
are stable in the range a — 4/a < A < —a + 4/a 
as was found for Stji- Nevertheless, there is a 
distinction between Sm and Sym fo r A < a with 
the dynamics of the YM field being different. 



(b) » YM ^ 

For the case with Y m 0, the YM field plays 
an important role in the dynamics of the universe. 
We find A = ±1 unless w = 0, for which we do 
not have any interesting dynamics. 

A = — 1 and 1 correspond to the case of the 
electric component dominance (B Y m) and that of 
the magnetic component dominance (B Y m), re- 
spectively. As we have already mentioned, the YM 
field always consists of both components. Hence 
these fixed points are reached only asymptotically, 
if they are stable. Just the same as the U(l) triplet 
fields, we find two fixed points for each case. How- 
ever, one of them with w = — 3AA is unstable. 
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Hence we discuss the other cases: 

(1) A = l (B YM ) 
4 



a + A ' 

_ a(a + A) - 4 
™ " (a + A) 2 
(2) A = -l (E YM ) 
4 

w = , 

a — A 

a(a - A) - 4 
(a - A) 2 



(4.27) 



(4.28) 



These hxed points correspond to the solutions with 
the magnetic component and the electric one found 
in §. IIIIA3l and lIIIA2| respectively. 

Without the non-linear interaction, these points 
were symmetric: they were related by the electro- 
magnetic duality and had the same stability prop- 
erties. The YM coupling skews the symmetry. 

For the case (1), the eigenvalues are given by 



8A 
a + A 
2(a - A) 



(4.29) 
(4.30) 



a + A 

and the two roots of the algebraic equation 

(a + X) 2 uj 2 + {a + X)(a + 3A)w 
+ [a(a + A) - 4] [4 + 3A(a + A)] = . (4.31) 

We find all eigenvalues are negative if and only if 
A > a and a (a + A) — 4 > are satisfied. Since 
this condition corresponds to the power-law infla- 
tionary solution with YM field, we can conclude 
that the power-law inflationary solution with mag- 
netic component dominance (Sym-I) is an attrac- 
tor. The difference from the U(l) multiplct case 
is that the solutions in the parameter range of 
< A < a, which are not inflationary, are no 
longer an attractor. We will discuss later which 
asymptotic state we find in this region. 

On the other hand, in the case with A = — 1, 
the eigenvalues are given by 



8A 
a — A 
2 



(4.32) 
(4.33) 



and the two roots of the algebraic equation 

(a - \) 2 uj 2 + (a- A) (a - 3A)w 
+ [a(a - A) - 4] [4 + 3A(a - A)] = . (4.34) 

We find that three eigenvalues are negative for the 
power-law inflationary solution as long as A < — a 



and a (a — A) — 4 > 0, but one eigenvalue wr, which 
corresponds to the perturbations of T (non-linear 
interaction term of the YM field), is always positive 
and does not depend on any parameters. This is 
the same behavior which we have seen in § 1111 B II 
As a result, this solution is unstable and the typical 
instability time scale is 0(1) e- folding time since 
the present time coordinate is A = ln(a/ao) and 
therefore T oc exp (u>rN). If the magnetic compo- 
nent is initially sufficiently small, we may find this 
power-law inflation solution by the electric com- 
ponents in the beginning, but the orbit leaves it 
just after 0(1) e-folding time. We conclude that 
the power-law inflationary solution with the elec- 
tric component dominance £ym _ I is unstable, con- 
trary to the Em . 



2. rv o 

Here we assume that Q Y m because V be- 
comes important when the YM field gives non- 
trivial contribution to the cosmic expansion. Note 
that A 2 = 1 (E YM or B YM ) is no longer a fixed 
point. We find the following two non-trivial fixed 
points, if A 2 > 6: 



w = w K 



(A± VA 2 - 16/3) 



The values at both fixed points can be described 
neatly by the deceleration parameter 



q 



(±) 



(4.35) 



as 



(±) 

YM 



o 
r (± 



( g (±))2(i + g (±); 



1 + ?(±) 



p 



(±) 



2(2 

1 



q 



(±) + i 



Note that £lv = at these fixed points and they 
are unique to the YM case. From (|4.35[) and the 
positivity of the density parameter £l YM > 0, we 
find 

l<q {±) <2. (4.36) 

Hence the power exponent of the scale factor p^> 
is 

1/3 <p {±) < 1/2, 

which is between a radiation dominant state and a 
stiff-matter dominant one. 
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The perturbative analysis is common to both of 
the fixed points A > y/6 (NA + ) and A < — a/6 
(NA-). We find the following eigenvalues; 

^r-A = g - 2 , 

which arc associated with the eigen vectors <5f2 Y M + 
(w/3)Sw and 5T/r~ [(q-2)(q+l)/8q]SA, respec- 
tively, and the two roots of the quadratic equation 

^ + (2 _ ?) , + to) =0 . 

I 



Since q is in the range of (|4.36p . if < A < a, 
we find all eigenvalues are negative, which means 
NA+ is stable. On the other hand, wn YM -ro turns 
positive when A < and we find NA_ is unstable. 
As a result, we find a stable fixed point in the 
parameter range of y/6 < A < a (NA + ), which 
partly takes care of the lost stability of Bym in 
the region A < a. 



We summarize our result for the SU(2) YM field 
in Fig. |3]and in Table HI 




FIG. 3: The parameter range for power-law solutions in the case with the YM field. The inflationary attractor 
solution with the magnetic component is found for a < A < — a + 4/a (Bym-I). On the other hand, the inflationary 
solution with the electric component, which is found in the range of a — 4/a < A < —a (_E Y m-I), is unstable. In 
the range of a — 4/a < A < — a + 4/a (SVm), we find the attractor solutions dominated by a scalar field. For 
inflation, we need an additional condition a < v2 (SVm-I). We also find new fixed points NA±, which exist only 
for non-Abelian gauge fields. 





SVm 






NA+ 


NA- 


existence 


a < V5 


a — 4/a > A 


A > —a + 4/a 


V6 < A 


A < -V& 


P 


2/a 2 


i(l-A/a) 


+ A/a) 


^tA + ^-ie/S) 


jt-(A- v/A 2 -16/3) 


SIym 





c(a-A)-4 
(<,-»■■> 


a(a + A)-4 
(«+A)2 


1(4 -A^ + A^A^- 16/3) 


|(4 -A^-A^/A 2 - 16/3) 


Qv 


1 - a 2 /6 


4-3A(a-A) 
3(c-A)^ 


4 + 3A(a + A) 
3(a + A)i! 








n K 


a 2 / 6 


8 

3(a-A) a 


8 

3(a+A)2 


1 (A- - 8/3 -AVA 2 - 16/3) 


|(A^ - 8/3 + A-^/A 2 - 16/3) 


stability 


a — 4/q<A<— a + 4/a 


always unstable 


a < A 


A < a 


always unstable 



TABLE II: The fixed points and their properties for the case with YM fields. Stability of SVm takes into 
account the fact that unstable V does not destroy the dominance of the scalar field. NA± cannot be inflationary. 
Inflationary conditions for the other points are the same as the U(l) triplet case. 



V. NUMERICAL STUDY 



From the above stability analysis, we find there are stable attractors if A > a (B YM A) or a > A > \/6 
(NA + ). We also find that a scalar field dominated universe (SVm), which is the same as the stable 
attractor in the model with a scalar field with an exponential potential (V = Vq exp(— a<p)), is stable in 
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the parameter range of a — A/a < A < —a + 4/a, even though the YM field does not necessarily settle 
down to its attractor state. As we will show here, it will oscillate in this scalar dominated background. 

We may also wonder what is the future asymptotic behavior for the other range of the coupling 
parameters a and A, i.e., A > a and A < Vo- Numerical calculations give us some insight into this 
question too. 

Numerical study also tells us strengths of the stable attractors. Since our stability analysis is based 
on the linear perturbations, we need numerical analysis to know how the attractor state is achieved from 
generic initial data. 



does not vanish. We indeed find the asymptotic 
value of the power exponent p is 3 instead of 2. We 



A. Numerical Analysis 

1. Stable attractors 

We begin with a small value of a for which the 
conventional power-law inflation is known to occur 
in the absence of gauge-kinetic coupling, namely 
a < V2. 

We choose the representative value to be a = 
1. We first performed the calculation for A = 2 
(see Fig. @|, which shows the conventional power- 
law inflation with an exponential potential (S VM -I 
). The YM field energy drops quickly. We find 
that the asymptotic power exponent of the scale 
factor is 2, which is consistent with the value of 
the conventional power- law inflation (p = 2 /a 2 ). 

1.0 p 



8 / 
h - 

4 I s 


p/5 












A/5 






V 


2 - 


10 
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FIG. 4: Inflation for a = 1, A = 2. This case obeys the 
usual cosmic-no-hair. 



When A > 3, our analysis suggests the power- 
law inflation assisted by the magnetic component 
of the YM field (B VM -V) is a stable attractor of 
the system. Figf5] confirms this fact as the density 
parameter for magnetic component stays constant 
(f2 Y M=constant and A = 1) while the scalar po- 
tential dominates the energy budget, which implies 
the universe undergoes accelerated expansion. An 
important difference between Figs[J] and [5] is that 
the acceleration is actually stronger when Dym 




FIG. 5: Occurrence of inflation assisted by the mag- 
netic field for a = 1, A = 5. Convergence to the inflat- 
ing attractor B YM T is clearly seen. 



deliberately chose the initial condition such that 
the scalar kinetic energy and the electric compo- 
nent are dominant over the others and the effect 
of YM coupling is significant. As shown in FigEl 
A approaches unity and T decays quickly whereby 
the system essentially reduces to the t/(l) triplet 
model. We find B VM -l asymptotically. 

Next, we take a negative A and confirm 
the electric-magnetic asymmetry for non-Abelian 
gauge fields. Fig|5] exhibits two different regimes. 
In the beginning, the electric energy density grows 
according to the linear instability caused by the 
strong gauge-kinetic coupling and the system is at- 
tracted towards the power-law inflation assisted by 
electric component of the YM field (E VM -T). 

During that period, however, T continues to 
increase and eventually destroys the inflationary 
regime at N ~ 10. The transient inflation E YM - 
I continues for 5^6 c-folding number, which is 
consistent with our evaluation given in §. 1111 B II 

After that, the universe is dominated by the 
scalar field while YM field is oscillating. In this 
case, since a is small enough to cause accelerated 
expansion by itself, this oscillation phase is also 
inflating. For comparison, we also show the plots 
with a smaller value of |A| (FigJTJ). The behavior is 
similar to Fig. [6j but there is no transient regime 
of E YM -l. When A is negative, from the instability 
of E YM -l, there is a peculiar behavior of rapid oscil- 
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FIG. 6: Inflation SVm-I with electromagnetic oscilla- 
tion via the electric YM inflation (Sym-I) for a = 
1, A = —5. The initial value of the ratio of the en- 
ergy density of the magnetic component to that of the 
electric one is 10 -8 . Between 5 < N < 10, 0,v is 
greater than its final value, which means the accelera- 
tion is stronger during that period thanks to the help 
by the electric component. As T gets to order unity, 
this regime is ruined and Qym decays while the fields 
oscillate. 
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FIG. 7: Inflation for a = 1, A = -2. While the dynam- 
ics of the universe is entirely dominated by the scalar 
field, gauge fields oscillate at late time. 



lation at late time between electric and magnetic 
components, which is not seen for positive A. 

Let us turn our attention to the supportive role 
of gauge fields in realizing inflation. We take a = 2 
for which inflation is impossible by the scalar field 
itself. With A = 5, we obtain Fig. [8] where A = 1 
in the future asymptotic state. The value of 
close to unity shows the expansion is accelerated, 
which can also be seen by the power exponent 
p > 1. As was investigated in the previous sec- 
tions, this is due to the interaction between the 
scalar and YM fields that transfers scalar field en- 
ergy to magnetic component of the YM field and 
slows down its rolling down the potential. 

Note that the velocity of the scalar field is given 
by </> = 2\nt/a, which is the same as the conven- 
tional power-law inflation. The difference is the 
values of total energy densities. The effective po- 
tential in the present model is given by the YM en- 




-0.2 J 



FIG. 8: Inflation assisted by the YM field (Bym-I) 
for a — 2, A = 5. Note that the asymptotic value of 
Qv is sufficiently large to cause accelerated expansion. 
Although the slope of the scalar potential is not flat 
enough to maintain the potential domination by itself, 
the magnetic component of the YM field also takes up 
the scalar field energy and helps realizing the inflation. 



ergy as well as the scalar potential V (Eq. (|3.6[1 ). 
which gives a larger Hubble expansion rate. As 
a result, the velocity with respect to the e-folding 
number N becomes slower as </>' = 4>/H. 




FIG. 9: Inflation for a = 2, A = — 5. Intermediately 
(2 < N < 6), the universe briefly inflates (i? Y M-I). 
Then F eventually dictates the dynamics of the YM 
field. 

For the U (1) gauge fields, the same type of infla- 
tion with non-flat potential E\j\-Y could have been 
seen for negative A because of the electro-magnetic 
duality. In the present non-Abelian case, however, 
a negative A drives not only electric component but 
also the normalized gauge coupling T, by which the 
inflationary regime is made transient and the final 
state contains mixture of electric, magnetic and 
scalar fields (FiglHJ) - During the transient phase 
of inflation supported by the electric component 
of the YM field (E VM -I), one can see the values of 
£V and f2 Y M being the same as the corresponding 
magnetic inflation. 

Finally, we confirm the stability of the new non- 
Abelian fixed point NA+ for A > (FigHUJ). The 
convergence is relatively slow and all the dynam- 
ical components undergo oscillations. In contrast 
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FIG. 10: Convergence to the non-Abelian attractor 
NA+ for q = 4, A = 3. It is distinct from the other 
plots in that T settles down to a constant value. 



to the other cases where T either diverges or dies 
away, this parameter region sees convergence to an 
attractor value, which necessarily means negligible 
Vly ■ Although it is not of interest in the context of 
inflation, it illustrates a distinct effect of the gauge 
coupling by forcing the potential term to vanish 
that would never happen in scalar-U(l) systems. 



2. Oscillation of the YM field 

The focus of this subsection is to understand 
the future asymptotic behavior of the system in 
the parameter region where the elementary fixed 
point analysis suggests there is no stable attractor 
solution. It turns out the nature of the dynam- 
ics in this regime is oscillation driven by the gauge 
coupling. 

Fig [TT1 shows the occurrence of scalar- YM oscil- 
lation as the future asymptotic state of the dynam- 
ical system for a = 4, A = —3. As is expected, the 
potential energy does not play a prominent role 
here, fly and f2 YM appear to converge to finite 
values although the numerical calculation has not 
been able to confirm it due to the computational 
difficulty caused by the rapid oscillation of A and 
the ever-growing T. 

The behavior is mostly the same for negative A 
regardless of A < —y/6 or not fFig |T2|) . The power 
exponent p of the scale factor is always slightly 
larger than 1/2. 

In contrast, A = s/E is a threshold value for pos- 
itive A since the YM fixed point becomes the at- 
tractor above it (see Fig |TU|) . Below the critical 
value, the asymptotic dynamics is rather analo- 
gous to the cases with negative A, but with a sig- 
nificantly smaller contribution of Qy- Convergence 
to an asymptotic value for fl YM can be seen more 
clearly here (Fig U3|) . The power exponent p of the 
scale factor is slightly smaller than 1/2 for A > 0. 




FIG. 11: Oscillation for a — 4, A = —3. As is expected, 
while r is smaller than unity, the system approaches 
the fixed point with the electric component Eym- Af- 
ter the effect of non-linear gauge coupling kicks in, the 
dynamics is irregular at the beginning. It appears the 
oscillation of f2y and Qym eventually die away, find- 
ing some asymptotic solution with the YM field oscil- 
lations. The power exponent p of the scale factor is 
slightly larger than 1/2. 



°\ T 1 




8 - '"\ / 




6 -/"\ jJ*w \ 






I ^ p 


4 " /~~^L 1 


1 n v 


2 - y \ 


MiMllilHiSSliill 


: a/5 2 / 

2 





FIG. 12: Oscillation for a — 4, A = —2. Qualitatively 
the same dynamics as A = —3. The only essential 
difference from Fig llll is the asymptotic value of Qv 
and D.YM- 



B. Asymptotic spacetime with the oscillation 
of the YM field 

From our numerical study, we find that the uni- 
verse still approaches some attractor spacetime 
but the YM field is oscillating for some parame- 
ter range, where we do not find stable attractors. 
In order to identify such an attractor by an ana- 
lytic approach, we assume that the time average 
of A, denoted by (A), does not change so quickly. 
We then discuss only three equations for w, f2 Y Mi 
and r, giving (A) = Ao (constant). From our nu- 
merical analysis, we find the following two typical 
asymptotic behaviors: 

(i) — > a finite value (A < 0) , 

(ii) fl v -> (A > 0) . 

r increases monotonically in our numerical study. 
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Appendix [5] Using this solution, we can take an 
average of A. However the background spacetime 
depends on Ao, which must be the same as the 
above averaged value (A) . Hence we need an iter- 
ative procedure to find the correct averaged value 
of Aq. In Fig. [Ml we present our result. 



Ao 



FIG. 13: Oscillation for a — 4, A = 2. One can see 
the exponentially decaying amplitude of oscillation for 
Sly M . In contrast to A < 0, the scalar potential contri- 
bution becomes completely negligible. 



We then do not consider the equation for T to find 
an approximate asymptotic solution. For the case 
(ii) , since there is the Hamiltonian constraint (|4.8|) , 
w and f2 Y M are not independent. 

We discuss the possible asymptotic solutions 
separately: 



1. Case (i) 

For the case (i), the dynamical equations are 

w' = - (6 - w 2 ) (a - w) + [2m- 3{a + AA )]f^ 
fl' YM = [-4 + w 2 + AkiAq + 4f2 YM ] f2 YM , 
where the reduced system gives a "fixed point" 
4 a(a + AA ) - 4 



a + AA 



(a + AA ) 



• (5-1) 



Using this "fixed point" , the equation for T is writ- 
ten as 



r' = 



- A 



AA 



r (> o) 



(5.2) 



This shows a monotonic increase of T, which is 
confirmed by our numerical calculation. 

The power exponent of the scale factor and the 
density parameter of the potential are given by 



AAn 



p = 



2a 



n v = 



3AA (q + AAp) 
3(a + AA ) 2 



.(5.3) 



From the positivity of density parameters, the 
following conditions must be imposed: 

a(a + AA ) > 4 , 3AA (a + AA ) + 4 > .(5.4) 

Once we know p and w of the background space- 
time, we can solve the YM equation as shown in 




FIG. 14: The averaged values of Ao (a), the power 
exponent p of the scale factor and density parameters 
1 (fi Y M! £V and Q,k) (b) for a = 4 and A < 0. 



As for the stability, we perturb the above two 
equations, whose eigenvalues are given by the two 
roots of the quadratic equation 



a + 3AA 
a + AA 
a(pt + AAo) - 
(a + XA y 



[4 + 3AA (a + AA )] =0. 



From the existence condition (|5.4[) , we find the 
following stability conditions: for a > 4/a/3, 



-a < AA < (AA ) ( -' 
where 



AA > (AAo)M ,(5.5) 



(AA ) 



(±) 



-a± v / a 2 - 16/3 



(5.6) 



while for a < 4/ a/3, we have only AAo > — a. The 
existence condition guarantees that two eigenval- 
ues are negative. As a result, this "fixed point" is 
always stable, although T diverges monotonically. 
We expect the universe in the parameter range of 
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(|5.4[) will evolve into this spacetime with the os- 
cillating YM field. Since AAo > 0, we find that 
p > 1/2, which is consistent with our numerical 
calculations. Note that for A = 0, by which we 
have a scalar field and Yang-Mills field without in- 
teraction, we find p = 1/2 as we expect. 

These approximate "fixed point" solutions seem 
to explain well our numerical results. 
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For the case (ii), using the relation £1 YM = 1 — 
^j 2 /6, we consider the following equation for w: 
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The asymptotic solution can be obtained as a 
"fixed point" in this sytem, which is 
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(5.7) 



It gives 



FIG. 15: The averaged values of Ao (a), the power 
exponent of the scale factor p and density parameters 
(Qym, £lv and SIk) (b) for a = 4 and A > 0. 



SIym — 1 77 A A 



o 1 



V 



4 + 3A 2 A2 



(< 1/2). 



(5.8) 
(5.9) 



In this background spacetime, we can also solve 
the YM equations as given in Appendix lAl Using 
this oscillating solution, we evaluate the averaged 
value Ao. However, since the background space- 
time depends on Ao, we have to find the correct 
value of Ao iteratively. In Fig. [T5l we show the 
result. 

Although the qualitative behavior coincides with 
our numerical result (for example, p < 1/2 and T 
increases exponentially.), it does not reproduce our 
numerical result quantitatively For instance, the 
asymptotic value of f2 Y M is ~ 0.7 in this approxi- 
mation, but the numerical value is ~ 0.4. A pos- 
sible source of discrepancy is that the oscillating 
time-scales for A and Q Y m are the same so that one 
cannot replace A by the constant averaged value 
Ao in the analysis of the dynamics of i7 Y M and 
tz7, even though the amplitudes of oscillations for 
those variables is dying away. 



VI. CONCLUDING REMARKS 

We have studied an SU(2) non-Abelian gauge 
field coupled exponentially to a scalar field with an 
exponential potential, while making a comparison 
with the U(l) multiplet case. 



We found that the power-law inflation with the 
magnetic component of the gauge field (£? YM T) is 
possible and it is an attractor of the present sys- 
tem, if A > a and A > — a + A/a. The transfer 
of scalar kinetic energy to the gauge fields through 
the gauge-kinetic coupling makes an inflationary 
solution possible even for a steep potential such as 
a > a/2, which is expected in the unified theories 
of fundamental interactions. 

On the other hand, the inflationary solution 
dominated by the electric component (i? YM -I) 
turned out to be unstable in contrast to the U(l) 
multiplet case. It can be a transient if the initial 
conditions are tuned. The attractor of the sys- 
tem is instead the conventional power-law inflation 
(S Y mT) if a < \/2. The YM field with a small am- 
plitude is oscillating in this background universe. 

We have also found new fixed points (NA±) in 
the parameter range of a/6 < A < a, which do 
not exist in the U(l) multiplet case. The fixed 
point NA + is an attractor, while NA_ is unsta- 
ble. We have also analyzed the non-inflationary 
regime, where the generic feature appears to be 
the oscillation of the YM fields (0±). 

We summarize our result in Fig. 1161 

One may wonder whether those isotropic in- 
flationary solutions are stable against anisotropic 
perturbations. Since there exist vector fields 
(A^), we usually find an anisotropic spacetime 
just as the case with a single U(l) gauge field. In 
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FIG. 16: The phase diagram in the parameter space 
of the present model. An inflationary phase is the at- 
tractor in the shaded regions. Besides, an attractor 
solution of the conventional sense exists for A > \/E 
(NA+). For the rest of the space, the nature of the 
dynamics is oscillation of the YM field. 



order to prove the predictive power of the scenario, 
we have to show that the FLRW universe is ob- 
tained as an attractor in anisotropic Bianchi cos- 
mologies. It is also interesting to know whether 
anisotropic inflation appears in a transient phase 
and its relic is observable or not. The study 
of Bianchi universe in the present model is in 
progress. 

Another important subject in the present model 
is a graceful exit from a stable inflationary uni- 
verse, including a reheating mechanism and a cal- 
culation of density fluctuations. In order to leave 
the power-law inflationary attractor that is a self- 
similar scaling solution, within the context of uni- 
fied theories of fundamental interactions, we may 
have the following possibilities, which may also 
work for the U(l) triplet case: 

(1) The moduli fixing: After a certain number of 
c-folds, if we can fix the moduli field <fi, the gauge- 
kinetic coupling vanishes. As a result, the inflation 
with magnetic component (£? YM -I) will end. 

(2) Hybrid-type Inflation: If Vo is not just a 
constant but depends on another scalar field a as 
Vo = ^-c , we find a dynamics approximated by 



the present scenario for the large value of a, and 
the end of inflation arrives when a gets small. 
(3) Decay of the VEV of YM field: The YM 
field may be coupled to other particles. Through 
such a coupling, the particles can be created quan- 
tum mechanically, which will reduce the YM vac- 
uum energy (/O Y m)01- The inflation assisted by 
YM field will eventually end. 

As for the reheating of the universe, for the cases 
(1) and (2), since we have a potential minimum 
around which a scalar field will oscillate, we then 
find the reheating of the universe. It is not clear 
whether we can find the hot Big Bang state via the 
particle production assumed in the case (3). 

The density fluctuations have been calculated 
for the case of U(l) triplet, which shows the lead- 
ing order effect of the background gauge fields is 
consistent with the current observational data [45$ . 
While YM field is expected to give qualitatively 
similar results at the linear order, there is an in- 
teresting prospect of generating non-Gaussianity 
through the famous chaotic behaviors that are pe- 
culiar to the non-Abelian gauge fields [H, Ugl, 
These subjects are to be investigated in future 
works. 
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Appendix A: Oscillation of the Yang-Mills 
field in the expanding universe 

In some numerical calculations, we have seen the 
YM field oscillates very rapidly while the back- 
ground spacetime evolves smoothly. If the energy 
density of the YM field is much smaller than that 
of the scalar field, the YM field does not contribute 
to the evolution of the universe. Even for the case 
that the YM field energy cannot be ignored, the 
oscillation of YM field may not directly affect the 
dynamics of the universe, but its mean value may 
contribute to the evolution of the universe. The 
different time-scales of the YM field oscillation and 
the evolution of FLRW universe may allow us to 
treat these two separately. Here we find such an 
oscillation of the YM field, assuming a given back- 
ground spacetime and evolution of the scalar field. 



Suppose the background is described by the fol- 
lowing power-law solution: 



a = a t p , 4> = wq-ZV + 0o , 



(Al) 



where p and zuq are constants, and N = ln(a/ao) is 
the e-folding time. The equation for the isotropic 
YM field in this background is given by 

A + P{1 + X ^A+ 2 4±^=Q. (A2) 
t al t 2 P y ' 

Changing the variables t and A to 77 and Z , which 
are defined by 



t 

A 



(?)' 



oo5ym fV 
. s 



V2 



— Xpmo s 1 3 



z, 



where 



3 - p(3 + Awo) ' 
we find the following equation for Z: 



(A3) 
(A4) 

(A5) 



d 2 Z 
dr] 2 



J 



[3(p-l) + 2\pw } 



2[3(p-l)+Apro ]s/3 



Z A 



(A6) 



Let us discuss the case where 77 increases as t in- 
creases, i.e., we assume that s > 0, or equivalently, 
3(p — 1) + XpniQ < 0. Hence this term in Eq. (|A6|) 
may drop as 77 — > 00 (t — > 00). Once we ignore 
the second linear term, we find a simple non-linear 
differential equation 



the YM field as 



d 2 Z 
'dq 2 



z 3 = . 



which solves as 



Z = Z cn [ Z ri; -j= 



(A7) 



(A8) 



where cn(x; k) is the Jacobi's elliptic function. 
Then the YM field is described in terms of the 
cosmic time t as 

A = aog 7 Z V A ^°/ 3 cn ( Zost 1 /*, -L) . (A9) 



V2 

Using this solution, we can evaluate the asymptotic 
behavior of the density parameter and the differ- 
ence between magnetic and electric components of 



0, 



x t 2-p(4+\m /3) 



A = 2cn 4 Z r] 



1 
71 



1 



(A10) 
(All) 



Since the above approximate solution contains the 
parameter s (|A5|) . which depends on the back- 
ground solution, there are the following two cases: 
(1) The background is controlled only by the scalar 
field. The YM field is oscillating in the back- 
ground, but its energy density is too small to affect 
the evolution of the universe. (2) The other case is 
that the averaged value of Ao as well as f2 Y M give 
an important contribution onto the background. 
In that case, we need an iterative procedure to 
find the correct averaged value Ao, as shown in 
the main body of the article. 

Here we present the averaged value of A and the 
properties of the asymptotic spacetime in the case 
(1). The results for the case (2) are given in §. 

EU 

For inflation driven by a scalar field, we have 
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p = 2/ a 2 and zuq = a. Then the condition 3(p — 
1) + \pvj < is 

A<A__J. (A12) 

This is always satisfied in the range we consider. 
The average of A must be taken in terms of the 
cosmic time t. We show our numerical result in 

Fig. EH 




FIG. 17: The averaged values of Ao for a = 1, v2 and 
2. It changes from —0.4 to 1 depending on the value 
of A. 



